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Abstract The mapping observations of CO J = 2 - 1, CO J = 3 - 2, 13 CO J = 2-1 
and 13 CO J = 3 - 2 lines in the direction of IRAS 22506+5944 have been made. The 
results show that the cores in the J = 2 — 1 transition lines have a similar morphology 
to those in the J = 3 — 2 transition lines. Bipolar molecular outflows are verified. The 
prior IRAS 22506+5944 observations indicated that two IRAS sources and three H 2 
masers were located close to the peak position of the core. One of the IRAS sources may 
be the driving source of the outflows. In addition, the H2O masers may occur in relatively 
warm environments. The parameters of the dense core and outflow, obtained by the LTE 
method, indicate that IRAS 22506+5944 is a high-mass star formation region. 
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1 INTRODUCTION 

Low-mass stars in molecular clouds are formed through the processes of collapse, accretion and out- 
flow (Shu et al. 1987). However, the processes for high-mass stars remain unclear because of the short 
evolution timescales and large distances of high-mass stars. Observational evidence suggests that young 
high-mass stars usually form in a cluster environment (Lada et al. 1993). Bipolar molecular outflows 
are commonly found around young high-mass stars (Shepherd & Church well 1996; Zhang et al. 2001). 
These outflows are generally much more massive and energetic than those from low-mass stars. Recent 
surveys with single-dish telescopes show that massive outflows are commonly associated with UC HII 
regions (Shepherd & Churchwell 1996; Zhang et al. 2001) and H2O masers (Elitzur et al. 1989). The 
H2O masers may originate in hot cores, and be excited by the shocks associated with outflows. However, 
how the outflows are driven is poorly understood. 

IRAS 22506+5944 has been proposed as a precursor of a UC HII region (Molinari et al. 1998). Two 
of the three H2O masers in IRAS 22506+5944 are found close to the outflow peak position (Jenness et al. 
1995). Wang (1997) imaged this source in the near-infrared J, H and K bands and found a star cluster 
within the core. It has magnitudes of 17.7, 14.5 and 11.7, respectively. Moreover, Zhang et al. (2001) 
detected bipolar outflows from an observation of CO J = 2 — 1 in this region. The bipolar molecular 
outflows are confirmed from observation of the CO J = 2 — 1 line (Wu et al. 2005). In addition, Su et 
al. (2004) reported interferometric observations in CO J = 1 - 0, 13 CO J = 1 - 0, C 18 J = 1 - 
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and the 3 mm continuum, as well as single-dish observations in CO J — 1 — and 13 CO J = 1 — 0. 
The IRAS 22506+5944 source was identified as the outflow driving source. 

In this paper, we report the first mapping observations of IRAS 22506+5944 in the CO J = 3 — 2, 
13 CO J = 2 - 1 and 13 CO J = 3 - 2 lines, as well as the mapping observation in the CO J = 2 - 1 
line. Due to the observed molecular lines at higher frequencies, we can attain higher angular resolution, 
which is critical for reducing the beam dilution and to identify the relatively compact core and outflows. 
Also, higher J transitions are relatively more sensitive to hot gases. Such hot gases are more likely to 
be physically associated with high-mass young stellar objects (Wu et al. 2005). 

2 OBSERVATIONS 

The mapping observations of IRAS 22506+5944 (R.A.(B1950) = 22 h 50 m 38.7 s , Dec(B1950) = 
59°44 / 58 // ) were made in the CO J = 2 - 1, CO J = 3 - 2, 13 CO J = 2 - 1 and 13 CO J = 3-2 lines 
using the KOSMA 3 m telescope at Gornergrat, Switzerland in April 2004. The half-power beam widths 
of the telescope at observing frequencies of 230.538 GHz, 345.789 GHz, 220.399 GHz and 330.588 GHz 
are 130", 80", 130" and 80", respectively. The pointing and tracking accuracy is better than 10". The 
DSB receiver noise temperature was about 120 K. The medium and variable resolution acousto-optical 
spectrometers have 1501 and 1601 channels, with total bandwidths of 248 MHz and 544 MHz, and 
equivalent velocity resolutions of 0.21 km s _1 and 0.29 km s _1 , respectively. The beam efficiency B e ^ 
is 0.68 at 230 GHz and 220 GHz. The beam efficiency £ eff is 0.72 at 330 GHz and 345 GHz. The for- 
ward efficiency F e ^ is 0.93. The mapping was done using on-the-fly mode with a 1' x 1' grid. The 
data were reduced using the CLASS (Continuum and Line Analysis Single-Disk Software) and GREG 
(Grenoble Graphic) software. The 80" resolution of t he CO J = 3 — 2 and 13 CO J = 3 — 2 data was 
convolved to 130" with an effective beam size of Vl30 2 — 80 2 = 102". The correction for the line 
intensities to the main beam temperature scale was made using the formula T m b = (F e ff / ^eff x ^a)- 

3 RESULTS AND ANALYSIS 
3.1 The Molecular Line Spectra 

Figure 1 gives the spectra of the different transition lines of CO isotopes at the IRAS 22506+5944 
position. Each spectrum shows broad line wings. The blue and red wings are asymmetrical and the 
spectral profiles are not Gaussian shaped. The observed parameters of the IRAS 22506+5944 source 
are summarized in Table 1 . In Table 1 , the larger full widths (FW) appear to indicate high- velocity gas 
motion in this region. The ranges of full widths are determined from the PV diagram. 



Table 1 Observational Parameters of the IRAS 22506+5944 
Source 



Name 






^mb 


FWHM 


FW 










(K) 


(kms- 1 ) 


(kms" 1 ) 


(kms" 1 ) 


CO J = 


2 


- 1 


12.8±0.1 


3.98±0.03 


11.11 


51.43±0.01 


CO ./ = 


3 


- 2 


11.9±0.1 


4.15±0.06 


13.63 


51.59+0.02 


13 CO J 




2-1 


4.1±0.1 


2.37±0.04 


6.69 


51.52±0.02 


13 CO J 




3-2 


4.0±0.3 


2.32±0.08 


7.22 


51.54+0.01 



3.2 Dense Molecular Core 

In Figure 2, all the structures of the molecular cloud core are presented and IRAS 22506+5944 has an 
isolated core. The cores in the J = 2 — 1 transition lines have a similar morphology to those of the 
J = 3 — 2 transition lines. Based on the results observed by other authors, there is a near-infrared source 
S4 which was detected by Wang (1997) as well as three H2O masers in this region. The IRAS sources 
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Fig. 1 Spectra at the central position of the region (0, 0) of the mapping 
observations. 



and the H2O masers are located closer to the core peak position traced by the transition J = 2 — 1, but 
a little deviation from that is also traced by the transition J = 3 — 2. The IRAS sources are associated 
with the H2O masers. This indicates that they are still undergoing activity in an early evolution stage 
(Felli et al. 1997). The H2O masers occur in relatively warm environments. Using the IRAS point- source 
catalog, infrared luminosity (Casoli et al. 1986) and dust temperature (Henning et al. 1990) are given 
by, respectively: 

L m = (20.653 x F12 + 7.358 x F 25 + 4.578 x F 60 + 1.762 x F100) x D 2 x 0.30, (1) 

_ 96 
TD = (3 + J) ln(100/60) - ln(F 60 /F 100 ) ' (2) 

where D is the distance from the sun in kpc; F12, F25, Fqq, and F100 are the infrared fluxes at four IRAS 
bands (12 /im, 25 /im, 60 fim and 100 /im), respectively. The emissivity index of the dust particles (/?) is 
assumed to be 2. The calculated results are presented in Table 2. 

The core parameters are calculated following the procedure of Lee et al. (1990). Assuming LTE, we 
write the column density as 

N = 2Ax 10 i 4 exp[/i£J(J + l)/fcT ex ] x (T ex + /iE/3K)xrxA^ cm _ 2 ^ 
J+l 1 — exp(— hv/kT ex ) 

where B is the rotational constant of the molecule, J is the rotational quantum number and v is the 
frequency of the spectral line. AV is a half full- width of speed, T ex is the excitation temperature and r 
is the optical depth. Here we assume A/'(H 2 )/iV(CO)=10 4 (Dickman 1978). 

The CO emission is always considered optically thick while the 13 CO emission is usually optically 
thin. Therefore, the excitation temperature T ex (CO) is given by (Garden et al. 1991): 



T ex (CO) = ^ 




^mb hv ( ( hv 



/ ' k V F VfcT b . 



(4) 
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Fig. 2 Integrated intensity maps of the core emission. In each map, the integrated range is 
from -53.26 to -49.68 km s -1 . The contour levels are 30%, 40%,..., 90% of the peak value. 
The red star represents the IRAS 22506+5944 source position, the green star represents the 
near-infrared source S4, and the filled triangles indicate the positions of the three H2O masers. 
The dot symbols mark the mapping points. 



Table 2 IRAS Fluxes and Derived 
Parameters 



Source 


Distance 


F12 


F25 


F 60 


i^ioo ^(^25/^12) 


lg(F 60 /Fi 2 ) T D 


^IR 




(kpc) 








(Ju) 


(K) 


(xlO 4 L ) 


IRAS 22506+5944 


5.1 


6.37 


34.74 


187.50 


295.10 0.74 


1.47 31.90 


1.38 



where Tb g =2.732 K is the temperature of the cosmic background radiation, and / is the beam filling 
factor. In addition, r( 13 CO) is given by (Garden et al.1991): 



r( i6 CO) = -ln<M 



hv \_exp(hv/kT ex ) — 1 exp(/iu/fcTb g ) — 1 

We also use another method to calculate t( 13 CO): 

T mb (CO) 1 - exp[-r(CQ)] 
T mb ( 13 CO) ~ 1 - exp[-r(i3CO)] ' 



(5) 



(6) 
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We assume the solar abundance ratio [CO]/[ 13 CO]=t(CO)/t( 13 CO)=89 (Lang 1980). The calculated 
results are listed in Table 3. From the table, we can see that the optical depths calculated by the above 
two methods are almost equal. Thus, we suggest that the abundance ratio [CO]/[ 13 CO] in this region 
is nearly the same as that in our solar system. Using the column density, the masses of outflow gas are 
obtained by: 

M = fiNn 2 S/ (2.0 x 10 33 ), (7) 

where the mean atomic weight of the gas /i is 1.36, and S is the size of the core region. The physical 
parameters of the core are summarized in Table 4. 



Table 3 Calculated Results of the Optical 
Depth 



Name 




Tl 


T2 


CO J = 2 


- 1 


33.55 


34.00 


CO J = 3 


- 2 


52.33 


53.59 


13 CO J 


2-1 


0.38 


0.38 


13 CO J = 


3-2 


0.59 


0.60 



The optical depth (ti) is derived based on the first 
method. (Y2) is derived based on the second method. 



Table 4 Physical Parameters of the 
Core 



Name 




r N(CO J = 2 — 1) 


Af(H 2 ) 


M 




(K) 


(xl0 17 cm- 2 ) 


(xl0 21 cm- 2 ) 


(xlO 3 M ) 


IRAS 22506+5944 


20.5 


33.55 7.07 


7.07 


2.03 



3.3 The Bipolar Outflows 

From Figure 3, bipolar outflows are clearly revealed in the IRAS 22506+5944 region. The red wing and 
blue wing lobes mostly overlap, which may be attributed to the axis of outflows being nearly parallel 
with the line of sight direction. The position-velocity (PV) diagrams in Figure 4, with a cut along the 
north-south direction, also clearly show high-velocity outflows. The outflows of the CO J = 3 — 2 are 
much clearer than those of the CO J = 2 — 1. Two IRAS sources are located close to the center of the 
outflows, so they may be the outflows driving the source, but we cannot distinguish which source is the 
driving source. 

Under the LTE assumption, the averaged column density of the outflows can be obtained by 
(Scovilleet al. 1986): 

N - 2 4 x 10 i4 exp[^J(J + l)//cT ex ] (T ex + ^/3K) f rdv 2 

N - 2AX1 ° JTi %xp(-W^ ex )y TmbX [l-exp(-r)] Cm * (8) 

We consider that the excitation temperature is uniform in the observed region. The excitation tempera- 
ture (T ex ) is 20.5 K from the calculated results of the cores. The r can be determined by Equation (6). 
The integral is carried out in the blue and red wing regions and the integrated ranges are determined 
from the PV diagram. 

In addition, we can derive the outflows' mass (M) from Equation (7). The momentum P and energy 
E are calculated by P = MV and E = MV 2 , where V is the mean velocity of the gas relative to the 
cloud's systemic velocity. A dynamic time scale can be determined by td = R/V, where R is the mean 
size of the outflows. The driving force is determined by F = P/td. The mechanical luminosity and the 
mass loss rate of the outflows are calculated using L mec h = E/td and M = Pj (tdV w ), where the final 
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Fig. 3 Left panel, outflow contours of the CO J = 2 — 1 are superimposed on the integrated 
intensity map of the core; the integrated range is from -57.07 to -53.26 km s _1 for the blue 
wing {black solid line) and from -49.68 to -45.96 km s _1 for the red wing (red dashed line). 
Right panel: outflow contours of the CO J = 3 — 2 are superimposed on the integrated 
intensity map of the core; the integrated range is from -58.96 to -53.26 km s _1 for the blue 
wing (black solid line) and from -49.68 to -45.33 km s _1 for the red wing (red dashed line). 
The contour levels for core and outflows are 30%, 40%, 90% of the peak values. 
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Fig. 4 Left panel: P-V diagram constructed from the CO J = 2 — 1 transition. Contour levels 
are 1, 2, 3, 4.5, 6, 6.5, 7.5, 9, 10.5, 12, 14,..., 25 K. Right panel: P-V diagram for CO J = 3-2. 
Contour levels are the same as in the left panel. The left and right vertical lines indicate the 
beginning of the blue and red wings, respectively. 



wind velocity is taken to be v w = 500 km s 1 (Marti et al. 1986). The calculated results are listed in 
Table 5. 



Table 5 Physical Parameters of the 
Outflows 



Name Wing 


N(H 2 ) 


M t d M 


F 




P 


E 


-^Mech 


IRAS 


(xlO 20 cnT 


- 2 )(M )(xlO 5 yr) (M© yr" 


x ) (MQkms- 1 


yr 


^(Mokms - 


^(xlO 46 ^ 


I) (Lq) 


22506+5944 Blue 


1.35 


18.7 4.14 8.0 x 10 


- 7 4.0 x 10 


-4 


164.6 


1.08 


0.22 


Red 


1.35 


18.6 4.38 8.2 x 10 


- 7 4.1 x 10 


-4 


179.6 


1.30 


0.24 
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4 DISCUSSION 

From Figure 1 and Table 1 , the full width of the CO J = 3 — 2 line is larger than that of the CO J = 2 — 1 
line and the line wings in CO J = 2 — 1 are broader than those in CO J = 1 — (Wu et al. 2005). These 
results suggest that the outflows may be caused by stellar wind sweeping up the surrounding materials 
in different intensity layers and the outflows traced by the higher J level CO lines arise from the warm 
layer closer to the central exciting star. 

Although the color indexes of the source satisfy the criteria of Wood & Church well (1989), they 
are not detected in the centimeter or millimeter continuum emission (Wu et al. 2005). Whether or not 
it is a UC HII region still needs to be determined by further higher sensitivity continuum observations. 
So far, the majority of authors have used constant values of r and T ex to estimate total mass. Because 
our observations were simultaneously made in CO and 13 CO, we can obtain relatively accurate values 
of r and T ex , as well as the total mass. From Table 4, the total mass of the core is 2.O3xlO 3 M ; 
the value indicates that IRAS 22506+5944 is a high-mass star formation region. Wang (1997) found 
clusters of stellar objects and an infrared jet in this region. Thus, we conclude that IRAS 22506+5944 
is associated with molecular cloud complexes. The IRAS 22506+5944 source appears to be a deeply 
embedded protostar. 

In addition, the contour maps and PV diagram of the IRAS 22506+5944 source clearly show bipolar 
outflows in this region. The total mass of the outflows is 37.3 M , and the total amount of momenta 
is 344.2 Mq km s _1 . Both parameters are much larger than the typical values in low-mass outflows. 
The larger outflow's mass suggests that the outflows' mass is not likely to originate from the stellar 
surface, but could be caused by the entrainment of ambient gas (Shepherd & Church well 1996). The 
IRAS 22506+5944 source is considered as the outflow driving source by Su et al. (2004), while Wu et 
al. (2005) consider S4 to be the outflow driving source. Two IRAS sources are located at the geometric 
position of the driving source of the outflow, and we also cannot identify which one is the driving 
source. We need to apply for observations with a much higher angular resolution telescope to confirm 
it. ^m/^Mech ^> 1, and Lir/CF < 1, where Z/Mech is the total mechanical luminosity, suggesting that 
despite the strong radiation from the central source, the radiation pressure still cannot supply enough 
force to drive such massive outflows (Wu et al. 1998). 

5 CONCLUSIONS 

We observed the CO J = 2 - 1, CO J = 3 - 2, 13 CO J = 2 - 1 and 13 CO J = 3 - 2 lines in 
the direction of IRAS 22506+5944. IRAS 22506+5944 has an isolated core. The infrared data indicate 
that the IRAS 22506 +5944 source appears to be a deeply embedded protostar. Bipolar outflows are 
further identified in this region. IRAS 22506+5944 or NIR source S4 may be the outflow driving source. 
The two IRAS sources are associated with H2O masers. The H2O masers occur in relatively warm 
environments. The parameters of the core and the outflows are derived; the derived values suggest that 
IRAS 22506 +5944 is a high-mass star formation region and the total mass of the outflows is larger than 
that in low-mass star formation regions. Compared with the extension of line wings in CO J = 2 — 1, 
CO J = 3 — 2 line wings extend further in velocity, suggesting that the outflows traced by the higher J 
level 12 CO lines arise from the warm layer closer to the central exciting star. 
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